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| We present a determination of the strange quark mass using lattice 

. QCD. Particular focus is put on the definition and renormalization of the 

mass. The latter is done non-perturbatively, using a recursive finite-size 
scaling technique. The hadronic regime of QCD, where the kaon mass is 
■ used as input of the calculation, is connected with the perturbative regime, 

where the strange quark mass can be translated into the MS scheme. A 
summary plot of the present lattice computations using dynamical (sea) 
O |' quarks is included. 
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^ ■ 1. Introduction 

h : 

Quantum chromo dynamics (QCD) as the theory of strong interactions 
between the Nf = 3 light quark flavours q = u,d,s (we treat the c, b, t quarks 
as infinitely heavy) is described through the Lagrange density 

I 

£QCD(go,m q ) = -— ^TriF^F^} + q(j ^(8^ + A^) + m q )q. (l.l) 
l % q 

The bare gauge coupling go and bare quark masses m q are the free pa- 
rameters of QCD. In order to make predictions, we need after a suitable 
regularization, e.g. on a Euclidean space-time lattice, to fix the free pa- 
rameters through a set of Nf + 1 physical quantities from experiment. For 
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Fig. 1. The hadron spectrum from quenched lattice QCD. 



example, we can take as experimental input hadronic quantities such as 



Experiment 
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and, using Eq. extract from the high-energy regime of QCD the so- 

called renormalization group invariant (RGI) parameters associated with the 
running renormalized gauge coupling and masses. In addition a number of 
other predictions can be made for phenomenologically relevant quantities. 
A crucial point becomes immediately clear: we need a tool to relate the 
hadronic, low-energy regime of QCD with the high-energy (> 2-10 GeV) 
regime, where perturbation theory in the gauge coupling applies. In the 
following we will demonstrate that such a tool is provided by the Euclidean 
space-time lattice. 

The predictive power of lattice QCD can be seen in Fig.^ which is taken 
from pQ. The data points represent various hadron masses m/, computed 
by the CP-PACS Collaboration The input used in the lattice simula- 
tions are m„, and m^, which fix the lattice spacing a, the mass of the 
degenerate up and down quarks, and the strange quark mass (a chiral ex- 
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trapolation is needed for the mass of the up and down quarks) . The hadron 
spectrum is extrapolated to the continuum limit a — » 0. Given that quark 
polarization effects have been neglected (the so-called quenched approxima- 
tion), the agreement between data points and experimental values marked 
by the horizontal lines is remarkable. 

In this article we will be concerned with the computation of the strange 
quark mass from lattice QCD. In the limit of vanishing quark masses (chiral 
limit) QCD possesses a large chiral symmetry. This symmetry is sponta- 
neously broken and the spectrum contains eight pseudo-scalar Goldstone 
bosons, which correspond to the observed eight lightest hadrons (7r's, K 7 s, 
rj). The latter are not exactly massless due to the non-zero quark masses. 
Quark masses are the explicit symmetry-breaking parameters 0] and they 
are treated as perturbations of the chiral limit in the framework of chiral per- 
turbation theory [3]. At lowest order the quark mass ratios m u jm^ = 0.56 
and rris/md = 20.1 can be determined from the pion and kaon masses [HlEj. 
Chiral perturbation theory, though, cannot determine the absolute scale of 
the quark masses; lattice computations are required for this. The strange 
quark mass can be determined almost directly through lattice computations; 
few assumptions are needed, which we will explain below. For the up and 
down quark masses, more difficult chiral extrapolations are needed. 

The value of the strange quark mass is required as an input for phe- 
nomenological predictions of the Standard Model, such as the CP-violating 
ratio e'/e. Also for physics beyond the Standard Model the quark masses 
are very important parameters. 

2. The strange quark mass from lattice QCD 

In the following we will use Wilson's formulation of lattice QCD, in- 
cluding Symanzik's 0(a) improvement. A review of principles of lattice 
computations can be found in Ref. [S]. The Wilson formulation explic- 
itly breaks the chiral symmetry through an 0(a) term in the Wilson-Dirac 
operator .9 . As a consequence the quark mass receives an additive renor- 
malization. However, it is well understood how the chiral Ward identities 
can be implemented on the lattice up to cut-off effects, which are 0(a 2 ) in 
the improved theory; for a short review see |10j . 

A renormalized quark mass on the lattice, translated into the continuum 
MS scheme, can be defined through 

m MS (/^) = Z m (g ,afj,)(m - m crit ) . (2.1) 

A definition that avoids the determination of the additive mass renormaliza- 
tion m cr it is through the partial conservation of the axial currents (PCAC) 
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relation, e.g. in the form 

M («7 M 75s)lat = ( m u + ""t s )latC"75s)lat + °( fl2 ) ' ( 2 - 2 ) 

where = («7 At 75s)i a t is an improved axial current ^Tj, P = (u'y^s)^ a 
pseudo-scalar density, and are the bare lattice current quark masses. A 
renormalized quark mass is defined from the PC AC relation Eq. ([2.2)1 using 
the renormalized axial current (^4r)^ = Z^A^ and pseudo-scalar density 
Pr = ZpP and comparing the bare and renormalized relations: 

(m u + m s ) m (fi) = Kim s ) kt Jl a(9o) . (2.3) 

The renormalization factor Z^{go) has been calculated non-perturbatively, 
using the chiral Ward identities [T^J. In order to compute Z P AS (go, o/f), we 
could apply perturbation theory. At one loop, 

2 P %,¥) = l + ||[(2/7r)ln(a^) + fc]+0( 5o 4 ), (2.4) 

where A; is a calculable constant that depends on the lattice regularization 
and the renormalization scheme employed. Bare perturbation theory is 
known to be unreliable at the coupling accessible to simulations go ~ 1 and 
the systematic error of the perturbative expansion is difficult to assess. A 
non-perturbative determination of Z-p avoids these problems. 

The running of the renormalized gauge coupling g and quark masses 
{fn q } depends on the renormalization scheme employed and is controlled by 
the renormalization group equations (RGEs) 



M^p = 0(g) , M^p 1 = r(g)m (l . ( 2.'> ) 



in terms of the (5 and r functions. They are in general non-perturbatively 
defined. Their perturbative expansions are 

P(g) ~ -f{bv + big 2 + b 2 f + ...}, (2.6) 

9-0 

r(g) ~ -fido + hf + ...}. (2.7) 

9-0 

If the renormalization conditions are imposed at zero quark mass ^3] the 
j3 and r functions do not depend on the mass. One such scheme, the MS 
scheme, is the one most commonly used for perturbative QCD computations 
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and there the (5 and r functions are known up to the 4-loop coefficients 63 
[TH ITS] and c?3 ^3 El • The RGI quark masses are defined through 



_ / _ 2 \-<V(2£>o) f /'SO) 

M g = m g (/x) (^26o5 (aOJ exp < - dz 



r(x) _ _tfn_ 
/3(x) box 



(2i 



Together with the A parameter, they are independent of the renormalization 
scale fj, and, in fact, any physical quantity, whose total dependence on fi 
vanishes, can be considered as a function of A and {M q }. The RGI masses 
are the same in all massless renormalization schemes. 

A straightforward implementation of the definition Eq. (|2.3f) on the lat- 
tice poses a scale problem. In order to compute Zp, the matrix element 
{0\P\K + ) is needed and for this the spatial lattice size L has to be large 
enough with respect to a typical scale like the kaon decay constant Fk, so 
as to avoid finite- volume effects. The renormalization scale \x has to be large 
enough to make contact with perturbation theory and at the same time \i 
has to be small with respect to the cut-off 1/a in order to avoid large cut-off 
effects. We get a set of inequalities ^H] 

11 11 
* F K 0.2GeV ^ /i lOGeV ^ ' 

which imply L/a > 50. But lattices that big cannot be simulated. One 
elegant solution to this problem is to take a finite-size effect as the physical 
observable, that is to identify L = with the only requirement L/a ^> 1. 
The renormalization scale [i is changed recursively in steps by factors of 2. 
Individually at each step the continuum limit can be taken. One therefore 
speaks of the recursive finite-size scaling (or step scaling) technique [T5j . 
We take the strange quark mass as an example: a running mass rn s (/j,) is 
defined 1 that runs with the system's size fi = 1/L. Starting at a reference 
maximal value L max , the system's size L is halved n times until contact with 
the perturbative regime is made; the RGI mass M s can then be extracted: 

£max = const. / F K = O(0. 5fm) : — ► L max m s (/i = 1/L max ) 

I 

FmaxTH'sifl' — 2/Z/ max ) 
I 

always a/j, = a/L <g 1 : 

I 

LragxfTl s (^fJ, — 2 /-t^max) 
perturbation theory j 

L ma ,x M s = # 



For a specific renormalization scheme implementing this idea, see below. 
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If the reference scale L max is expressed in units of the kaon decay constant 
Fk, the result is a number for M s /Fx, which is scheme-independent. 

A particularly suitable renormalization scheme to implement the finite- 
size scaling technique is the Schrodinger functional (SF). It has been defined 
for QCD in Refs. (201111] and is the "work horse" of the ALPHA Collabora- 
tion. QCD is formulated in a box where Dirichlet boundary conditions are 
imposed on the fields at Euclidean time xq = and xq = T. They provide 
an infrared cut-off proportional to 1/T to the frequency spectrum of quarks 
and gluons, making it possible to perform simulations, and hence impose 
renormalization conditions, at zero quark mass. A renormalized gauge cou- 
pling g(L), running with the scale /j, = 1/L, is defined through a variation 
of the effective action with respect to a change of the boundary gluon fields 
|22l I23| . Keeping g{L) fixed, while changing the bare lattice parameters, 
is equivalent to keeping the physical size L (and hence the renormalization 
scale /i) fixed. 

The recursive finite-size scaling technique in the SF has been applied 
to compute the running of the renormalized coupling with Nf = 2 flavours 
of massless quarks |23j . The non-perturbative j3 function is shown by the 
points in Fig. [21 together with a non-perturbative fit. It deviates from the 
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Fig. 3. The running quark mass in the SF scheme for Nf — 2 flavours of massless 
quarks. Plot taken from |27| . 



3-loop (5 function for ase{L) = 5 2 (L)/(4vr) > 0.25. 

The renormalization factor Z-p can be defined and computed in the SF 
scheme |24 (l25| l26]. The running of the renormalized quark mass is extracted 
from it through 



m 



■SF 



(M) 



m SF 



W2) 



lim 



Z 



SF/ 



,2L/a) 



Zf(g ,L/a) 



(2.9) 



H=l/L 



which is independent of the quark flavour q. The results for Nf = 2 flavours 
of massless quarks have been presented in detail in Ref. [27] and are shown 
by the points in Fig. |21 The running starts at the low energy n = 1/L max ~ 
0.4 GeV corresponding to the leftmost point and extends to high energies, 
where perturbation theory can be safely applied to extract the RGI mass. 
The curves are obtained by using the perturbative expressions for the t / (3 
functions at the loop order indicated in the legend. In this case perturbation 
theory works surprisingly well down to small energies. For the running 
gauge coupling, instead, at asF(-^max) = 0.367, clear deviations between 
non-perturbative and perturbative curves can be seen in Fig. |2 

In the following we summarize the computation of the strange quark 



8 



PROC PRINTED ON FEBRUARY 2, 2008 



mass done in Ref. [21] ■ The steps involved are visible in the equation 

-MS/ n /™W\ M ( M ^ F i ^ ^A(go) /oim 
mf /x 2 = — — x — - — x (m s ) ut — r . 2.10) 

The scale in the MS scheme is conventionally ^2 = 2 GeV. The first factor 
on the right-hand side of Eq. (12.10(1 is the perturbative (PT) evolution in 
the MS scheme. The second factor is the non-perturbative (NP) running in 
the SF scheme that starts at the hadronic scale /xi = 1/L max . The third 
factor is the renormalized (at the scale /ii) strange quark mass in the SF 
scheme and it has to be computed at several lattice spacings to take the 
continuum limit. This part involves in principle 

1. simulations of u, d and s quark flavours at their physical masses, and 

2. setting the overall scale of the simulations, i.e. determining the lattice 
spacing a in MeV -1 . 

Because of computer cost and technical difficulties, simulations of two 
light and non-degenerate u and d quarks with an additional s quark have not 
been possible so far. There are very promising developments reported in 28 , 
which show that these difficulties can be overcome in the near future. For 
the time being, the strategy of [23] is adopted. Simulations are performed 
with two degenerate flavours at a reference mass m re f. The latter is tuned 
so that the pseudo-scalar (PS), made of two identical flavours, has a mass 

mps(m re f,m re f) = rn K = - (jn 2 K + + m |-o) QCD = 495 MeV . (2.11) 

The subscript "QCD" means that an estimate of the electromagnetic effects 
has been subtracted from the experimental numbers [20] in order to obtain 
a pure QCD kaon mass, as we have on the lattice. 

The lattice spacing could be set by computing on the lattice the kaon 
decay constant, a number aFx, and dividing it by the experimental value of 
Fk, which can be obtained from the decay rate of K + — > Instead, 
in the present computation, the lattice spacing is set through the scale ro 
extracted from the static quark potential [20]. Data for r^/a as well as 
pseudo-scalar masses needed to tune m re f in Eq. (|2.11|) are taken from 31 , 
where they are available at three lattice spacings in the approximate range 
0.092-0.071 fm. The data for r^/a are extrapolated to zero quark mass |23| . 
The phenomenological value is ro = 0.5 fm, obtained from potential models. 

The results for the reference RGI quark mass show a strong dependence 
on the lattice spacing 27 o . No systematic continuum extrapolation is there- 
fore performed, instead a continuum estimate is given: 

M ref = 72(3)(13)MeV, (2.12) 
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where the value and the first error are the result at the smallest lattice 
spacing, the second error is a systematic one and is the difference with 
the value at the largest lattice spacing. The latter dominates. The result 
Eq. ()2.12|) is consistent with the one obtained in the quenched (iVf = 0) 
theory 2 |29j . It is assumed to hold in the Nf = 3 theory as well. In order 
to translate Eq. (|2.12|) into a value for the strange quark mass, the Gell- 
Mann-Oakes-Renner formula is used 

m\ = (M + M s ) Brgi = 2M ref B RG i, (2.13) 

where M = \{M U + Mj) and -Brgi is a constant of the chiral Lagrangian. 
Together with M s /M = 24.4(1.5) 32 , Eq. (t2~T3l can be solved for M s , and 
using 4-loop running in the MS scheme yields 

mpGu = 2GeV) = 97(22) MeV . (2.14) 

The m^ s values for the three lattice spacings are displayed in Fig. 0] (filled 
red squares). 

3. Conclusions 

In Fig. 0] we show a summary of the lattice determinations of the strange 
quark mass with dynamical (sea) quarks at finite values of the lattice spacing 
a. This is an update of the plot of Refs. [2Z1E3]- There are points obtained 
with JV f = 2 dynamical quarks HZlEIlElESlES]andAr f = 2 + l [S3 EH1 
In the legend we indicate the gauge action/fermion action employed and 
whether non-perturbative or perturbative renormalization has been used. 
The dictionary for the gauge actions is: W: Wilson action; I: Iwasaki action; 
LW: 1-loop tadpole improved Liischer-Weisz action. The dictionary for the 
fermion actions is: W: Wilson action; CW: Wilson-clover action; KS: Asqtad 
staggered action. Except the Wilson fermion action used in jHH], all the other 
lattice formulations have 0(a 2 ) discretization errors. The non-perturbative 
renormalization is done in the SF or the RI-MOM 32j scheme. The physical 
kaon mass is always used as input to fix the strange quark mass, 

There is a systematic effect due to the use of perturbative renormal- 
ization, which leads to values of the quark mass smaller than those for 
non-perturbative renormalization. Figure 0] also shows the impact in the 
quark mass values of |38| when the perturbative renormalization factor Zp 
in Eq. (|2.4|) is evaluated at two loops |4U| : the mass values increase by 
15%. There is good agreement between data from different non-perturbative 

2 Also for the A parameter there is no significant difference between Nf = and Nf = 2 

m 
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■ ALPHA, N f =2, W/CW, nonpert. renorm. (SF) 

► ALPHA, Nf =2, W/CW, nonpert. renorm. (SF), Z£ 

• CP-PACS, N f = 2, l/CW, pert, renorm. 
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T JLQCD, N f = 2, W/CW, pert, renorm. 

A QCDSF/UKQCD, iV/ = 2, W/CW, nonpert. renorm. (RI-MOM) 

♦ SPQcdR, N f = 2, W/W, nonpert. renorm. (RI-MOM) 

O CP-PACS/JLQCD, N f =3, l/CW, pert, renorm., stat. err. only 

£ HPQCD/MILC/UKQCD, N f = 3, LW/KS, pert, renorm. 1-loop 

'~ □ HPQCD/MILC/UKQCD, JV> = 3, LW/KS, pert, renorm. 2-loop 

0.005 0.01 0.015 0.02 0.025 0.03 



Fig. 4. Summary plot of the strange quark mass from dynamical simulations with 
N { = 2 and iV f = 3. 

renormalization procedures and, in fact, between them and the Nf = re- 
sult |29| marked by the dotted lines. We also mention the most recent 
compilation of the strange quark mass from QCD sum rules jljj, which 
quotes m s {2 GeV) = 99(28) MeV, in agreement with the range of the lattice 
results; see also [321 H31 H2] . 

In conclusion, Fig. 0] demonstrates that non-perturbative renormaliza- 
tion is essential for a reliable lattice determination of the strange quark 
mass and that data at smaller lattice spacing(s) are needed for a systematic 
continuum extrapolation. 
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